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ff^ ' Abstract 



We construct a mass function for the Schwarzschild-de Sitter black hole spacetime everywhere in 
between the black hole and the cosmological event horizon. The mass function is positive definite, 
continuous and monotonically increasing with the radial distance from the black hole horizon and 
can be related to the geodesic motion. Finally, choosing a suitable timelike hypersurface in between 
the two horizons, we study particle creation for a massless quantum scalar field in this spacetime. 
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1 Introduction 

In general relativity, an important concept is that of a mass function. The principal features of any 
mass function should be threefold. First, it must be defined with respect to a timelike translational 
Killing vector field; second, one should be able to relate the mass to the geodesic motion for a 
Newtonian interpretation and finally, the mass function must be a positive definite quantity. It is 
the third criterion that makes it a difficult problem to define a mass function because one cannot 
define a satisfactory notion of conserved gravitational energy-momentum tensor unless one goes 
to the asymptotic region (see e.g. [I] and references therein). Only an approximate notion of 
gravitational Hamiltonian density can be defined perturbatively and locally but the positivity of 
this quantity is far from obvious. 

For asymptotically fiat spacetimes a gravitational mass can be defined in several ways. The 
simplest is the Komar mass. This is proportional to the surface integral of the derivative of the 
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norm of the timelike Killing field and thus is related directly to geodesic motion. In general the 
Komar integral will be positive definite only if the matter energy-momentum tensor Tab satisfies 
the strong energy condition (SEC) : (Tab - \Tgab) ^"C'' > 0, for any timelike C° 1]. 

The second is the Arnowitt-Deser-Misner (ADM) O |3l |4] formalism. In this approach a gravi- 
tational Hamiltonian density is defined with respect to the background timelike Killing vector field 
in the asymptotic region and the integral of this Hamiltonian density is computed in the asymptotic 
region. This integral is interpreted as the gravitational mass. 

It is known from the Raychaudhuri equation that geodesies would converge in a mass distribution 
only if the latter satisfies the SEC [Tl[5l[6]. Also, it is known that the SEC implies the weak energy 
condition, i.e. the positivity of the energy density. Using these two facts a third approach to 
define gravitational mass and to prove its positivity was developed in [7l IB] for asymptotically flat 
spacetimes. 

The positivity conjecture of the ADM mass was first proved in [HI HU]. Soon afterward, a re- 
markable proof of the positivity of the ADM mass appeared in [11] using a spacelike spinor field 
on a spacelike non-singular Cauchy surface. This result was generalized for black holes in asymp- 
totically flat or anti-de Sitter spacetimes in [T2]. The A < spacetimes usually have well deflned 
asymptotic structure or inflnities which are accessible to the geodesic observers. The references 
mentioned above consider explicit asymptotic structures of such spacetimes at spacelike infinities 
which are uniquely Minkowskian or anti-dc Sitter. In fact the positivity of the ADM mass for 
A < physically reasonable spacetimes admitting spin and well defined asymptotic structures is 
well understood so far. 

But recent observations suggest that there is a strong possibility that our universe is endowed 
with a small but positive cosmological constant A [131I14]. We note that, since a positive A violates 
SEC, it repels geodesies (see e.g. [121 HH HZl UHl [H]) and hence the first and the third of the 
methods mentioned above to define gravitational mass do not seem to be applicable in this case. 
Also, the known exact stationary solutions with A > (see e.g. [20]) usually exhibits an outer 
horizon, namely the cosmological event horizon. The tiny observed value of A (of the order of 

10~^^m~^) sets the length scale of the horizon to be O which is of course very large but 

finite. If a black hole is present, it will be located inside the cosmological horizon and the spacetime 
is known as de Sitter black hole spacetime. The cosmological event horizon acts in these spacetimes 
as an outer causal boundary [21] , beyond which the timelike Killing vector field becomes spacelike 
and communication is not possible along a future directed causal path thereby ruling out any precise 
notion of asymptotics. 

The very first approach to define mass in de Sitter black hole spacetimes appeared in |21| . 
where a mass function was defined on the black hole and cosmological horizons using the integral 
of their respective surface gravities. The variation of this mass function gave a Smarr formula. 
In [22] ■ metric perturbation was considered in a region far away from the black hole but inside 
the cosmological event horizon where the background spacetime is de Sitter. A local gravitational 
energy momentum tensor was constructed and with respect to the background de Sitter timelike 
Killing field the mass of the perturbation was defined. This perturbative approach has much 
similarity with the usual Hamiltonian formulation of general relativity. We shall adopt this approach 
explicitly in this paper. For asymptotically Schwarzschild-de Sitter spacetimes the mass in this 
asymptotic region with respect to the de Sitter background was found to be M, i.e. the mass 
parameter of the Schwarzschild-de Sitter metric. The spinorial proof of positivity of ADM mass 
for asymptotically flat spacetimes [11] was generalized later in [23] |2l] to show that the mass thus 
defined in the sense of [22] with respect to the background de Sitter spacetime is indeed a positive 
definite quantity. 

What do we wish to do then? Since there exists no preferred asymptotic region for defining 
mass we shall define this everywhere in the region between the black hole and the cosmological 
horizon, keeping in mind that firstly it must be a continuous positive definite quantity. Secondly, 
since positive A corresponds to a positive energy density, the A part of this mass function should 
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increase monotonically with radial distance from the black hole horizon. Thirdly this mass function 
must be related to the geodesic motion in order to have a satisfactory Newtonian interpretation. 

We outline the basic scheme now. We shall set c = /cb = G = ?i = 1 throughout. The metric 
signature will be taken mostly positive (—, +, +, +). If not otherwise mentioned, repeated indices 
will be summed over. Let us consider the Schwarzschild-de Sitter spacetime written in spherical 
polar coordinates. 



2M 
r 



Ar 



2M Ar 



(1) 



For 3M\/A < 1, this spacetime has three Killing horizons. 



ru = 
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cos 



-cos-^ (smVa) + - 
3 V / 3 



rc 



Va 
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rv = - (rn + rc) 



rn is the black hole event horizon and rc > th is the cosmological horizon, whereas ry is unphysical. 
For 3Ma/A < 1, it is easy to see that 2M and rc ^J^■ 

As mentioned earlier, the gravitational mass of the perturbation over the de Sitter background 
has been defined and computed earlier in [52] in a region where ^1 — ^ ^ and it turned out 
to be M. Instead, we divide the region between the black hole and the cosmological event horizon 
(ru < r < rc) into three consecutive regions of perturbation 

f2M Ar2\ , . / Ar2\ 2M . . . / 2M\ Ar^ 
+ ^ Region I , 1 ^ > (Region II , 1 > — 



1 » f — + —1 (Region I), M - _ 1 » _ (Region II), (l - — 1 » — (Region III) 



where in the first inequality each term on the right hand side is much smaller than unity, the term on 
the right hand side of the second inequality is much smaller than each of the terms on the left hand 
side and so for the third inequality. In this sense the three regions are distinct. These three regions 
can respectively be interpreted as perturbations over background Schwarzschild, Minkowski and de 
Sitter spacetimes. For A ^ 10~^^m~^ and M ranging between the extremes 10'*m to lO^m, it can be 
easily checked that the above regions exist and are merged smoothly in between. We shall compute 
the gravitational mass for each region following [22\. Then by defining the gravitational mass of 
the background in a suitable way, we shall construct the total mass function of the spacetime. The 

perturbation scheme obviously fails near the horizons, ^1 — _ -> 0, and we shall treat 

these regions separately. 

After we define the mass function, we will connect this to timelike geodesic motion. We shall 
compute its variation to get a Smarr formula. We shall further present a calculation to determine 
the spectrum of particles created in this spacetime for a massless quantum scalar field in a suitable 
timelike hypersurface in a region between the two horizons. 

The paper is organized as follows. In the next section we shall construct the total mass function 
and connect this to geodesies. In Section 3, we shall discuss the Smarr formula and the particle 
creation. Finally we discuss our results. 

2 Construction of the mass function 

We start by considering the A- vacuum Einstein equations 

Rab - ^Rgab + Agab = 0, (4) 



where Rab and R are respectively the Ricci tensor and scalar computed from the metric. Let us 
assume we can find a region in between the black hole event horizon and the cosmological horizon 
where the metric gab can be decomposed in a background g^^J and a perturbation jab over it 

gab = gfb + lab, (5) 

where |7Qf,| ^ Idl^bl- ^^'^ basic scheme described in 122 can be outlined as follows : define a local 
'gravitational energy-momentum tensor' T^f^ which consists of quadratic and higher order terms 
of the 7's, whereas the Einstein tensor consists of g^^J and terms linear in jab- Let vi'^'' denotes the 
metric compatible covariant derivative on the background g^J ■ Then the 'energy current' T^^''^'''-*" 
is conserved with respect to the background, i.e. v['''' (^T^^-'a^'"-'''^ ~ 0, where ^C^)" is the timelike 
Killing field corresponding to the background 

V(0)^(0) + vW^) « 0. (6) 

Then one computes the flux of the energy current over a closed spacelike hypersurface S. We shall 
apply this method to compute the gravitational mass of the perturbations in different regions of 
Eq. ©. 

Let us start by considering region I of Eq. i.e. linear perturbation over the Minkowski 
spacetime 

Itt = ( — + — )'>'■= I ~ + ~ ) ' 7ee = = 7^^, (7) 
and the components of the background Killing vector field, 

0, 0, 0}, 0, 0, 0}. (8) 

The Ricci tensor and scalar reads 

Rac = + \ [V("'^ (V(°)7ee + ^fhae - ^fhac) - Vi") {^^"^^cf + ^^^h " ^^"^'jc, 

+0{j') + ..., 



0(7') + ..., (9) 



where the trace is defined with respect to g'ff^ 17 — jabg^^'^'^'' ■ Also, for the Minkowski background 
which is a A = vacuum, the Einstein equations become identities 

r':j - Ir^"'9^^ = 0. (10) 

Now we use Eq.s Q and PH]) to expand Einstein's equations dH) with — 0. We shift the 
0(7^) and other higher order terms to the right hand side of Eq. Q which define the gravitational 

energy-momentum tensor (rj^'j^^ — Kgab 



V(°"^V(")7,. + V(°)'^Vr 7.. - V(")'^v(")7a. - (V(°)''V(°)^7..) ffi",^] = ^^T^^' - ^9ab, (H) 



where jab = lab - ll9ab ■ 

The gravitational mass {Mq) is defined as the integral of the 'energy current' {T^ — Kgtb) 
over a spacelike hypersurface S orthogonal to the timelike Killing vector field. 



(G) 



Kgtb) e'rfS, (12) 



4 



where 'i' corresponds to the direction of the timehke Kilhng field. Following (22], we obtain the 
following expression for the energy current from Eq. after a little algebra, 



1 



tb 



^ 167r 



y(0)d 



tbcd 



(0)c 



HtbcM'^'e 



where 



^^abcd — I i/c i< . V . i( . 



7ca 76d - 5cd Taf, - 9ab led + 5bd lea 



(13) 



(14) 



Habcd is antisymmetric under the interchange of (a, d) and (6, c). Then since in Eq. (jl3p the 
indices are fixed a = t = b, i.e. timelike, the indices (d, c) must be spacelike. Therefore we can 
convert the integral in Eq. ([T2|) into a surface integral 



Ma - 



IGtt 



(15) 



where 'dS' denotes the volume element of a closed 2-surface, i.e. the boundary of E of the region 
of interest E. Since Schwarzschild-de Sitter spacetime is spherically symmetric, the closed surface 
is a 2-sphere. Also, Eq. ([7]) gives 



7 = labg 



{Q)ab 



0. 



(16) 



We now explicitly evaluate Eq. (IT5|) using Eq.s (O, ([5]) and (ITB)) . The covariant derivative on the 
background Killing vector field is V(o)d^b ~ dd^b ~ T^j^^e where we keep in mind that d and b are 
summed over as in the equation. Since Habcd is antisymmetric under the interchange of a and d, and 
a is timelike, d must be spacelike. Keeping in mind 5 = i, it is clear that V(o)d^& is non-vanishing 
only when d = r. We also have c = r in the second integral of Eq. ([15]). But Eq.s (jH]), ([7]) and 
give Httrr — 0. Thus the second integrand in Eq. (|15p is identically vanishing. Now expanding the 
covariant derivative in the first integral we find that the only non- vanishing term is —g'^'^^lelfd- 
Since the direction 'd' is along r, we finally get 



Ar3 

Mg = M+ — 
6 



(17) 



Next we consider perturbation over the de Sitter background, i.e. region II : (^1 — ^ in 

Eq. (JSj. This was explicitly done in [22]. The perturbation and the background Killing field are 
the following 



and 



(0) 

9tt 
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Ar2 



2M 
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2M 



Ar 
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c(o)'^^{i, 0, 0, 0}, e;.' 



(0) 



Ar2 
~3~ 



0, 0, 



(18) 



(19) 



The calculation of the mass of the perturbation is essentially the same as before. The only difference 
we have to remember is that the de Sitter background we are considering now is a A-vacuum 



Kb - 2^ 9ab + ^9ab - LI- 



(20) 
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The calculation of the mass of the perturbation, using Eq.s ^ and ([20]) and following the 
method described above gives 

Mg = M. (21) 

Thus the mass of the perturbation with respect to the background de Sitter spacetime is given by 
the mass parameter of the Schwarzschild-de Sitter spacetime. 

Finally we consider region III in Eq. ([3]), i.e. perturbation of the Schwarzschild background by 
a A term. The background and the perturbation are 

3 3(1-2M) 



0, 0, 0}, ^^(l- — ), 0, 0, OL (23) 



and 



We follow the same procedure described after Eq. (IT6l) . The second integral in Eq. (fT5|) can be 
shown to be vanishing as earlier and evaluating the first integral we get 

Ar^ 

Mg = -r- ■ (24) 
6 

The mass functions of Eq.s ([T7| . (|2T|) . (|24|) are clearly positive definite, but are not continuous. In 
order to achieve a continuous mass function we introduce the notion of the total mass function in 
the following way. First we note that for the Minkowski background, the background curvature is 
identically vanishing. Therefore we take the background mass to be zero. For the Schwarzschild 
background we define the background mass to be the Komar mass, 

MB^-^<f^aii°^dS-, (25) 

where the integration is performed over a 2-sphere. For Schwarzschild spacetime, Mb = M any- 
where inside that particular perturbation region. 

For the de Sitter background we treat the —Agf^ term appearing in the Einstein equations 
as the energy- momentum tensor (SttT^^) corresponding to the cosmological constant. Thus the 
corresponding energy current becomes 

T.\e'''' = -^9'fJ(''''. (26) 
an 

But we have from the unperturbed A- vacuum equation ()20p . 

Ta\^'''' = -^Rllik^'^' = ^V(")'^v(°)e\ (27) 

using the Killing identity. Thus the gravitational mass Mb corresponding to the background de 
Sitter vacuum is given by 

Mb = J T^tS^^^'dS = i- y v(°)'^V;,°)e^°^dS. (28) 
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Since ^(o)"^ is a timelike coordinate Killing field, the index d is spacelike above, as can be seen by 
antisymmetrizing the covariant derivative on So we can convert the integral in Eq. (j28p into 
a surface integral over a 2-sphere to get 

MB^^/Vi")^")^^'^^^. (29) 

We now combine the gravitational mass of each region with the respective background mass to get 
a total mass function 

U{r, M) = Mb + Mq = M + ^ , (30) 



which is positive definite and monotonically increases with r, thereby encompassing the satisfaction 
of weak energy condition by positive A. 

The perturbation theory fails as one moves close to the horizons 

2A/ Ar^X 

l----)^0, (31) 

since we cannot naturally decompose the spacetime in two parts among which one contains a back- 
ground timelike Killing vector field. We can bypass this difficulty as the following for 3Mv^ <C 1, 
which is of course the physically reasonable situation. Let us first consider the region infinitesimally 
close to the cosmological event horizon. Instead of taking Schwarzschild-de Sitter spacetime, we 
assume a solution of the kind 

2M Ar^\ ,o ( 2M Ar^" 



ds2 = -(l he — jdt'^ + il he — I dr^ + r^dO"^ + sin^ 



(32) 



where e is a dimensionless infinitesimal positive constant. It can be checked that this metric satisfies 
Einstein's equation. The effect of e is to make the cosmological event horizon larger and the black 
hole smaller. Now, near the cosmological horizon the quantity is itself an infinitesimal. We 

choose e ^ O ( — ) . It is then obvious that we can choose 
/2M \ 2M-re 



7tt = e , 7rr = — -—2,706 = = 700, (33) 
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and the background timelike Killing field same as described in Eq. (|19p . With this, we repeat 
the same procedure described earlier and we take e — )■ at the end of the calculation to get Eq. (1501) . 
Similar analysis can be performed infinitesimally close to the black hole event horizon (with 

e ~ C (^""P^^ ) to get the same mass function. Thus we have constructed a mass function everywhere 
in the region between the black hole event horizon and the cosmological horizon. The mass function 
is positive definite, continuous and monotonically increasing with r, i.e. the more we move away, 
the more positive A-energy we get inside. 

How is this function C/(r, M) = + really a mass? More precisely, is it meaningful to 

call this the mass inside the sphere of radius r? One way to answer this question is by studying 
timelike geodesies in the Newtonian approximation for ([Ij . For any geodesic in this spacetime, four 
conserved quantities i?, Li, L27 -^3 can be defined respectively to the timelike Killing field and the 
three generators of a 2-sphere. The first is the conserved energy and the rest can be regarded as the 
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components of the orbital angular momentum. Then the timelike geodesic motion can be mapped 



to an effective 1-dimensional non-relativistic motion of a unit mass test particle with energy ^i?^ 
as in the Schwarzschild spacetime |1], 



if^+V^r, = (34) 

where the effective potential ipir, L) is given by 

, , ^, 1 / 2M Ar2\ /L2 \ 1 U{r,M)L'' U{r,M) 



where 



Ll + Ll + Ll = (e'^ + sin^ 6*) , (36) 



and [/(r, M) is the mass function given by Eq. ([50]). For a Newtonian interpretation we note that 
for sufficiently large r, we have 



1 AL2 uir^M) , £^ 



The above equation clearly indicates that a geodesic motion far away from the black hole in the 
Schwarzschild-de Sitter background can be thought of as a motion under a positive definite, position 
dependent mass function U (r, M) , along with the usual repulsive centrifugal barrier term and some 
constants. 

Alternatively, we can write down the geodesic equation explicitly 

(p■x^' „ dx'' dxP , , 

where r is the timelike parameter along the geodesic, and take its non-relativistic limit. We set 
H — r and consider region I of Eq. ([3]) and take the limit t —> t. The equation becomes 

i^.(f-^)-^j^i^-..(-^)-a;-^ 

using Eq. The above equation clearly shows again the justification of calling U{r,M) a mass 

function. 



3 Smarr formula and particle creation 

We shall now compute the variation of this local and total mass function (|5D|) subject to the change 
of the black hole mass parameter M and keeping A fixed. For A = this gives rise to the idea of 
black hole thermodynamics |25j [26l [27l [28l [29] (see also [30] for a vast review). 
The area of the black hole horizon (rn) is given by 



^H=47rr2, (40) 



which, using Eq. 1^ we rewrite as 
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Now we write the mass function U{r, M) in terms of two new variables: the black hole horizon area 



Ah and volume V 
the mass function, 



enclosed by an arbitrary sphere of radius r on which we have defined 



UiAK,V) 
the variation of which gives 

6UiAa,V) = 



4A [AB_y fA^ 
' 3 V1671- / VIBtt 



2A 



(^h)^ 



(167r)2 2(167r)2 (Ah)' 

The surface gravity kh of the black hole horizon is given by 



Ay 

"87 



Stt 



KH = 



M Arn 



combining which with Eq.s PO)) . P5)) gives 

5U{An,V)^^5A^ + -^5V. 
We can do a similar calculation with the cosmological horizon in Eq. PO)) to write 

Stt Stt 



(42) 



(43) 



(44) 



(45) 



(46) 



where Aq and kc are the cosmological horizon's area and surface gravity respectively. The term 
^ appearing in either of Eq.s (|45|) . (|46| can be interpreted as a negative isotropic pressure due to 



A 

8^ 



A 



UaUb 



{UaUb + gab) 



(47) 



where Ua is a unit timelike vector along the worldline of fluid particles and the energy density p 
and the isotropic pressure P of the fluid have been identified with ^ and — ^ respectively in an 
orthonormal frame. Combining Eq.s (j45p and (I46p for the same volume V we get a Smarr formula 
involving horizon parameters only 



khMh + KcSAc = 0, 



(48) 

which was derived earlier in \21\ I31j , which we have rederived using our local mass function. Note 
that since both kh and kq are positive, the Smarr formula shows that when the area of the black 
horizon increases, the area of the cosmological horizon decreases and vice versa, an immediate 
consequence of the expressions of the horizon lengths, Eq.s ([2]), which state that for 6M > 0, th 
increases and rc decreases and vice versa. It is clear then unlike asymptotically flat or anti-de 
Sitter spacetimes, one cannot interpret ^ as the entropy of the spacetime. Instead one defines a 
'total' entropy S — ^ri+Ac j3Yj Por further discussions on this we refer our readers to [32l [33l |34] 
which consider the variation of A as well. 

Since kh kc for rn 7^ Eq.s (HSI) . show that unlike asymptotically flat or anti-de 
Sitter spacetimes, there can be no unique thermodynamic interpretation in the Schwarzschild-de 
Sitter spacetime unless one separates the two horizons using a thermally opaque membrane. In this 
situation one expects two different thermal equilibrium states of temperatures ^ , ^ corresponding 
to the black hole and the cosmological horizons respectively. On the other hand, the variation of 
the 'total' entropy combined with the variation of the total mass function U gives 



5U 



27r (kh + Kc) 



Ac 



A 



-SV, 



(49) 
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derived earlier in [32 using the global mass function defined on the horizons. Thus with respect to 
the total entropy one might expect an 'effective equilibrium' temperature T = 2Tr(K^+K.c) ' 

In the following we shall derive the different thermal equilibrium states using Kruskal patches 
for an eternal Schwarzschild-de Sitter spacetime corresponding to each of the horizons (i.e. corre- 
sponding to Eq.s (|45|), (|46l) ). 

It was shown in 35 using canonical quantization and Bogoliubov transformations that a stellar 
object undergoing gravitational collapse in an asymptotically flat spacetime to form a black hole 
creates Planckian distribution particles at late times. For massless quantum fields this distribution 
can be measured at the future null infinity, and found to have a temperature where kjj is the 
surface gravity of the black hole future horizon. Later this result was rederived using path integral 
quantization |36| . This most remarkable result, known as Hawking radiation, was further justified 
by the renormalization of the quantum energy-momentum tensors (see |371 138] and references 
therein). We further refer our reader to |39l|4^ for excellent reviews on this subject. 

For an eternal horizon, there is no scenario for gravitational collapse and there exists a past 
horizon, in addition to the future horizon. It was shown in [21| using path integrals that for 
an eternal Schwarzschild-de Sitter spacetime, the black hole and the cosmological horizon create 



thermal particles with temperatures ^ and ^ respectively. In 



particle creation on both 



the horizons was studied, showing that there can be non-thermal spectra. In [41' particle creation 
by a Schwarzschild black hole sitting within a Friedmann-Robertson- Walker (FRW) universe was 
studied. The FRW universe could be a global de Sitter space itself. See also [42] for study of 
particle creation in Schwarzschild-de Sitter spacetime via complex path analysis. 

We shall derive in the following the particle spectra within both the horizons, rjj < r < rc, 
using canonical quantization. We consider a massless minimally coupled scalar field ip moving in the 
Schwarzschild-de Sitter spacetime, and ignore any backreaction. Employing the usual separation 
of variables, 'ip{t, r, 6, <j)) — e~'"* ^"''"^''^^ Yimid, (p), the equation of motion for a single mode becomes 

d^firn{r,t) d\f,,Jr,t) 2M 



where we have abbreviated f{r,t) 
dr 



- 1- 



■In 



7(r 



Ar2 
" ~3~ 
and r 



M 



/im(r,t)-0, (50) 



is the tortoise coordinate defined by. 



2kh 















-^In 




+ ^ln 


^ + 1 




2kc 


rc 


2ku 





(51) 



where ku = dr [I — ^ 1- ) . Thus — ^ =Foo as r — > th, rc respectively. Since we are 

\ / r— ru 

ignoring backreaction, we have I <^ M . Clearly, to find and work with the exact mode solutions 
of Eq. ((50)) seems a formidable task and hence we shall have to choose a suitable region to work 
in. But before we go into that, we digress briefly to construct suitable coordinate systems for the 
Schwarzschild-de Sitter spacetime. There are two coordinate singularities located at th and rc, 
therefore we require two Kruskal-like patches to remove them. We define the usual outgoing and 
incoming null coordinates (u, v) as 

u = t — ri,, i) = t + r^. (52) 

By writing the metric ([T]) in terms of u and r, it is easy to find that u — > ±oo as r — !■ ru, tq 
respectively along an incoming null geodesic, whereas by writing it in terms of v and r gives v — > =Fcx) 
as — > th, rc respectively along an outgoing null geodesic. In terms of the null coordinates (u, v) 
the metric becomes 
2M 



r 

1 



r 

1 



^ + 1 



dudv 



(53) 



where r as a function of (u, v) is understood and can be found from Eq. 
Kruskal null coordinates for the black hole event horizon as 

1 _ _ 1 



. We now define the 



(54) 
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so that u — > 0, — oo as r — !■ 7'H I respectively, and v 0, + oo as r — >■ rjj, respectively. The 
Kruskal null coordinates for the cosmological event horizon can be defined as 



1 

— ( 



(55) 



so that 



-oo, as r -> th, rc respectively, and — > — oo, as r — > rn, fc respectively. To 



summarize, the ranges of the various null coordinate systems are 

— oo < u < oo, — oo < w < oo, — oo<u<0, 0<U< oo, < u' < oo, 



oo < TJ' < 0. (56) 



Clearly, there will be both outgoing and incoming mode solutions for the field equation. Since 
{u, u') and (tJ, TJ') are respectively functions of (u, v) only, we have modes in terms of these null 
coordinates. 



Ipin 



, t t , -t-t I -'t-'t 

biV, + blvj = hvi + blvj = 6-iJ- + fefwf , 



(57) 



where {ui, Vi), {ui, Vi) and (it-, «■) are modes corresponding to the coordinates (u, f), (u, v) 
and (m', v') respectively. The index 'i' corresponds to all discrete and continuous indices. The 
complex quantities etc. are expansion coefficients and as in flat spacetime, they are interpreted 
as creation and annihilation operators associated with respective modes. Note that for an eternal 
Schwarzschild-de Sitter spacetime, there should be only ingoing modes at r = rn and only outgoing 
modes at r = rc. In that sense the mode expansions in Eq.s ()57p . which contain both ingoing and 
outgoing modes, are complete. 

The creation and annihilation operators are defined to satisfy the commutation relations 



flj, a] 



[a,:, a,- 



= (5y , [ci, aj] = = 



al, a] 



6,,, bl 



bi, b. 



b,] ^ = 



= Sij, bj] = = 



a-, a'J 



-'t -'t 



b^, bj 



S,j, [a-, a'j] = 
The inner product of the modes (w^, Vi) are defined as 

{ui, Uj) = \ {^{VaUj) - Uj(Vaul)) dS'' = (5y , (Wi, Vj) 

{u„ u]) = 0= (wj, w|). 



b[, b'j 



bl b]_ 

I't T't 



(58) 



(«|(Va«,)-«,(V„4)) dS'' 



where S is suitable hypersurface and the direction 'a' is along its normal. In an asymptotically flat 
spacetime, one chooses E to be the past null infinity. But as we discussed earlier, in presence of 
a de Sitter horizon, infinities are not accessible to an observer located within that horizon. So we 
have to choose S differently here. 

In order to do that let us look at Eq. ((50| . The quantity — ^) can be written as -^dr {l — - 
Since the norm (l - ^ - ^) vanishes at two points ^rc, it must reach at least one maximum 
in between. This is given by = (^^)^, which ranges between lO^^m to lO'^'^m for rn ^ 10'*ni to 



lO^^m, and for the observed A 



A 

1-52^-2 



10 ''•^m ^. In any case, ^ 



1 



2M 



3 



1 near this maximum 



and we are therefore in region I of Eq. ([3]). Also, since by our assumption / ^ M, Eq.s (|50)) . (l52t . 
(j54p , (j55p show that around this maximum we have usual outgoing and incoming positive frequency 
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plane wave modes, 



u{u!,l,m)^ Yim{0, (/)), v{uj,l,m)r^ YimiO, (/>), 

u'{ujj,m)^^= Yi^{9, q^), T;'(w,/,m)^-= 0), (60) 

along with their negative frequency counterparts. Let us now define the Bogoliubov transformation 
coefficients (see e.g. [38] for details) and consider the outgoing u and u modes first, 

(Ui, Uj)=a,j, (ui, u]^ =/3ij- (61) 



Let us now consider the equality between the first two mode expansions in the first of Eq.s ff57\ 
We use Eq.s ([59l) . ([61]) and the commutation relations ([58| to get 



a^kal^ - l^ikPlj = <5y , a^kPlj - ^Ija^k = 0. (62) 

Subject to these relations and the commutations, one can then take the inverse transformations 

a.i = aijttj - ^Jjo]. (63) 

If |0) denotes the vacuum associated with the {u, v) modes, then the Kruskal observer will 'see' 
particles in |0) in the i-th mode, 

(0|a|a^|0) = |/3yf (no sum on i). (64) 

Note that we are doing this in a region where particles and number operator can have well defined 
meaning. With all the above equipments, our sole task is now thus to compute l3ij. We note that 
on any r = constant hypersurface, dt = d{t — r^,{r)) = du = e'^^'^du, using Eq.s ([54]) . There will 
be a dr^ coming from the normal direction of the hypersurface volume element. But dr^e^^'^^ — 
iwe"*"^" = — 9„e~*'^" = ^g-KH" ^g_^-ium(u)'^ ^ which means for these modes dtdr^ = dud^. Putting 
these in all together it is straightforward to calculate 



ik 



du, (65) 



where all the constants including those arising from summation of the discrete indices and angular 
integral have been dumped into the constant k. We are yet to choose the limit of the above 
integration. An observer sitting on that hypersurface will see null waves moving to and fro. Since 
on or in the vicinity of that surface the solution has plane wave solutions, the observer will see the 
phases of those waves to be constant. Since t varies on that surface, we exhaust all the possible 
values of those constant phases. Thus we choose the limit of u in Eq. ((65|) to be — oo to (Eq. ([SS)) ). 
With this, the integral in Eq. (I55|) looks exactly the same as in asymptotically flat spacetime and 
clearly we have been able to computationally map our problem to the asymptotically flat spacetime 
case [35l [39l [40] , and bypassed all the possible difficulties of a de Sitter spacetime regarding non- 
trivial mode solutions, asymptotics or even redshift effects. Following [351132], we now evaluate the 
integral. Analytically continuing this to the complex plane, and treating It = as a branch cut one 
obtains 

\au:,u>'\'^ ^ e~ \l3uj,uj'\'^ ■ (66) 
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Then from Eq. (|62p we get 



J duj'5{LJ^u')=-k' j dw'(l-e^) (67) 

where k' is some positive constant arising from summation of the discrete indices and angular 
integrations. We note that for kh — > 0, /3 — > and there is no particle creation. 

Thus the Kruskal observer will 'see' the |0) vacuum is filled with thermal distribution of outgoing 
particles around r — 

- , (68) 

e "H — 1 

with temperature = f^- Note that unlike the collapsing scenario we do not have any greybody 
effect here, because we are merely comparing two outgoing modes. Similar analysis for incoming 
modes {v, v) shows that there are incoming thermal particles with the same temperature too, 
showing that there will be thermal equilibrium. This justifies the thermodynamic interpretation of 
Eq. (I45p . The ASV term can be interpreted as the infinitesimal variation of r around r = r^. 

Similar calculations with incoming (w, v') modes show that the p rimed Kruskal observers will 
detect thermal radiation at temperature Tq = This justifies Eq. (|46| also as a thermodynamic 
equation. We also note that if we set kq = 0, there will be no particle creation for the cosmological 
horizon. 

For fermionic field the commutations are replaced with anticommutations and Eq.s (|62p is 
modified with a '+' in place of ' — '. This will give Fermi-Dirac distribution with the same respective 
temperatures. 

Let us summarize the results now. We have constructed a positive definite, local and continuous 
total mass function everywhere in between the two horizons of the Schwarzschild-de Sitter black 
hole spacetime. The A part of the mass function increases monotonically with radial distance, 
which is in correspondence to the positivity of the vacuum energy density. We have also related 
the mass function to the geodesic motion. We have rederived all the thermodynamic relations by 
varying this mass function. Finally we have computed particle creation in a suitable region of this 
spacetime to verify those thermodynamic relations. 
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